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ABSTRACT 

We report new Spitzer Space Telescope observations from the IRAC and MIPS 
instruments of the young (~ 3 Myr) a Orionis cluster. The populous nature of 
this cluster makes it a good target for statistically-significant studies of disk 
emission as a function of mass. We identify 336 stars as members of the cluster 
using optical and near-infrared color magnitude diagrams. Using the spectral 
energy distribution (SED) slopes in the IRAC spectral range, we place objects in 
several classes: non-excess stars, stars with optically thick disks (like classical T 
Tauri stars), class I (protostellar) candidates, and stars with "evolved disks"; the 
last exhibit smaller IRAC excesses than optically thick disk systems. In general, 
this classification agrees with the location expected in IRAC-MIPS color-color 
diagrams for these objects. We find that the evolved disk systems are mostly 
a combination of objects with optically thick but non-flared disks, suggesting 
grain growth and/or settling, and transition disks, systems in which the inner 
disk is partially or fully cleared of small dust. In all, we identify 7 transition 
disk candidates and 3 possible debris disk systems. There appears to be a spatial 
extension of infrared excess sources to the north-east, which may be associated 
with the young (< 1 Myr) embedded cluster NGC 2024. As in other young stellar 
populations, the fraction of disks depends on the stellar mass, ranging from ~10% 
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for stars in the Herbig Ae/Be mass range (>2 M ) to ~35% in the T Tauri mass 
range (1-0.1 M ). We find that the disk fraction does not decrease significantly 
toward the brown dwarf candidates (<0.1 M ). The IRAC infrared excesses 
found in stellar clusters and associations with and without central high mass stars 
are similar, suggesting that external photoevaporation is not very important in 
many clusters. Finally, we find no correlation between the X-ray luminosity and 
the disk infrared excess, suggesting that the X-rays are not strongly affected by 
disk accretion. 

Subject headings: infrared: stars: formation — stars: pre-main sequence - 
open cluster and associations: individual ( a Orionis cluster) — protoplanetary 
systems: protoplanetary disk 



1. Introduction 



Observations have shown and theory predicts that disks around young stars are an 
inescapable consequence of sta r formation, given angular momentum conservation during 



proto stellar cloud core collapse (jHartmannll2005al : I Ward- Thompson et al.ll2005l : lAndre et al. 



20001 ) . As the resulting protoplanetary disks age, their excesses at near-infrared (NIR) and 
mid- infrared wavelengths decrease. Previous studies have indicated th at the timescale for 



this evolution in disk e mission is strongly dependent on the stellar m ass (ILada &: Ladalll995 



Muzerolle et all 12003 : ICalvet et al.l 12004 : ISicilia-Aguilar et all l2005f ) ranging from 5-7 Myr 
for objects in the mass range of the T Tauri stars ( TTS, w ith spectral type K5 or later; 
Haisch et al1l200ll : lHartmannll2005bl : iHillenbrand et al.l 120061 ) to <3 Myr for objects in the 



mass range of the Herbig Ae/Be (HAeBe) stars (spectral type F5 or earlier: iHernandez et al. 
2005h . 



The observed decreases in infrared emission can result from: grain growth to sizes much 
larger than the wavelength of observation; dust settling to the disk mi dplane, which reduces 
the " flaring" of the disk and thus the amount of ener gy radiated (IKenyon fc Hartmann 



19871 : iDullemond &: Dominikl|2005l : iD'Alessio et al.ll2006l ): clearing of small dust particles by 
large bodies in the disk; or some combination of these processes. S ince grain growth and 
settling to the mid-plane oc cur fastest in the (warmer) inner disk (jWeidenschillingl Il997l : 
Dullemond &: Dominikl 120041 ). the disk emission is expected to decline faste r at short wave- 
lengths than at long wavelengths, and the re is observational evidence for this (ISicilia-Aguilar et al 
20061 : lHartmannll2005bl : lLada et allhood ). 



Of particular interest are "transition disk" objects, which have an inner, optically thin 
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disk region, possibly produced by planet clearing, combined with an outer, optically thick 
disk JCalvet et all I2OO2I I2OO5I : b'Alessio etaD l2005al : iMegeath et all l2005ah . Prior to the 
complete or nearly-complete clearing seen in transition disks, one would expect to observe 
optically-thick but reduced infrared emission from the inn er disk, for which there is some 
evidence (e.g., ISicilia-Aguilar et al.ll2006l ; lLada et al.ll2006l ). However, additional large sam- 
ples of stars are desirable to make the result more statistically si gnificant, as stars of the 
same mass and age show a wide range of disk emission properties ( iFurlan et al.ll2006l ). 



The cr Orionis cluster is of particular interest because it is reasonably near and relatively 
populous, making possible statistically-significant studies of disk properties as a function of 
stellar mass. Moreover, a Orionis has an age ~ 3 Myr at which one might expect the 
beginnings of disk evolution t o become evident. It is part of the Orion OBlb sub-association 



which has an age of 1.7-7 Myr ljWarren fc Hesserlll978l ; iBrown et al.lll994l . Il998l ; iBriceno et al. 



Walter et al. 



2005). Low-mass members of the a Orionis cluster were first reported by 
( 119971 ). who found over 80 X-ray sources and spectroscopically identified more than 100 low 
mass, pre-main sequence (PMS) members lying within 1 degree from the star a Ori. The 
distance calculated by Hipparcos for this star (352+g5 6 ) agrees, within the uncertainties of 



Hipparcos, with the Hipparcos distance calculate d for the overall population of the stellar 
subassociation OBlb (439±33; IBrown et al.l 119981). which is statistical l y more reliable. The 



estimated age of t he cluster is 2-4 Myr (IZapatero Osorio et al.l 120021 ; lOliveira et al.l 12002 



Sherry et al.l 12004] ; iFranciosini et al.ll2006l ). Since the a Orionis cluster is relatively near and 



the reddening tow ard the center of the cluster is low (E(B-V)=0.05 mag; IBrown et al.lll994 



Bejar et al.l Il999l ). this stellar cluster is an excellent laboratory to study young stars in a 



obiects. like brown dwarfs and free-floating planets (e.g.. Barrado v Navascues et al. 2001, 


2003; 


Beiar et al. '. 


999. 


2001. 2004: Burningham et al. 2005: Caballero 2005: Kenvon et al. 


2005: 


Oliveira et a 


. 2006: Scholz k Eisloffell 2004: Sherrv et al. 2004: ZaDatero Osorio et al. 


2002; 


Walter et al. 


1998 


)• 



The unprecedented sensitivity and spatial resolution provided by the Spitzer Space Tele- 
scope in the near- and mid-infrared windows are powerful tools to expand significantly our 
understanding of star and planet formation processes. In this contribution, we analyze the 
near- and mid-infrared properties of stars in the a Orionis cluster, ranging in mass from 
HAeBe stars to the substellar limit. 

This paper is organized as follows. In §2] we present the observational data. In S|3] we 
describe the selection of cluster members. We analyze the observations and describe the 
results on and give our conclusions in $5j 
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Observations 



2.1. Infrared photometry 

We have observed a field of ~1730 arcmin 2 on the o Orionis cluster usi ng the four 



chann els (3.6, 4.5, 5.8, and 8.0 \xm) of the InfraRed Array Camera (IRAC, iFazio et al. 



20041 ); 90% of this field was also observed using; th e 24 /zm band of the Multiband Imaging 
Spectrometer for Spitzer (MIPS, iRieke et al.ll2004l ). 



The IRAC observations presented here were taken on October 9, 2004. The field of 
view was covered by a 9 x 10 position mosaic (20" overlap), with three dithered exposures 
at each position. Images were obtained in High Dynamic Range (HDR) mode, whereby a 
short integration (1 second) is immediately followed by a long integration (26.8 seconds). 
Standard Basic Calibrated Data (BCD) products from version S13.2 of the Spitzer Science 
Center's IRAC pipeline were used to make the final mosaics. Post-BCD data treatment 
was performed using custom IDL software (Gutermuth et al. 2004, Gutermuth et al. 2006) 
that includes modules for detection and correction of bright source artifacts, detection and 
removal of cosmic ray hits, construction of the long and short exposure HDR mosaics, and 
the merger of those mosaics to yield the final science images. 

MIPS observations were made using the medium scan mode with half-array cross scan 
overlap, resulting in a total effective exposure time of 80 seconds. The 24 /zm images were 
processed using the MIPS instrument team Data Analysis Tool (DAT), which calibrates the 
data and app lies a distortion cor rection to each individual exposure before combining into a 
final mosaic ([Gordon et al.ll2005l ). 



Figures [T] and [2] show the IRAC color and the MIPS 24 /im images of the a Orionis clus- 



ter. T he lower-right panel in Figure [2] shows a map of dust infrared emission (ISchlegel et al. 



19981 ) illustrating the location of the cluster in the OB lb sub-associ ation, which can be 
traced by the ringlike structure centering at RA~84 and DEC~-2 ( jBriceno et al.l 12005 



Hernandez et al.ll2006l ). In this panel, the boxes show the positions of the IRAC (small box) 
and MIPS (large box) scan regions. The isocontours are estimators of gala ctic extinction 
(for A y=0.2, 0.4, 0.6 & 0.8 mag) from the 100 /im map of dust emission (see, ISchlegel et al. 



19981 ). The galactic extinction and the background at 8 and 24 /im is quite low compared 
with other regions with similar ages, which improves the precision and sensitivity of the 
photometry. 

Point source detections were carried out individually on each IRAC band using Phot Vis 
(version 1.09), an IDL GUI-based photometry visualization tool developed by R. Gutermuth 
using the DAOPHOT modules ported to IDL as part of the IDL Astronomy User Library 
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( [Landsman! Il993l ). More than 10,000 sources were detected in at least one Spitzer band. 
We extracted the photometry of these objects using the apphot package in IRAF, with an 
aperture radius of 3". 7 and a background annulus from 3.7 to 8". 6. We adopted zero-point 
magnitudes for the standard aperture radius (12") and background annulus (12-22". 4) of 
19.665, 18.928, 16.847 and 17.391 in the [3.6], [4.5], [5.8] and [8.0] bands, respe ctively. Aper 



ture corrections were made using the values described in IRAC Data Handbook (jReach et al. 



20061 ). A preliminary list of 1682 objects (hereafter Sample 1) was created by selecting those 



objects having photometric measurements in all IRAC bands with errors less than 0.15 mag- 
nitudes. We obtained point source photometry at 24 /im with IRAF/ daophot point spread 
function fitting, using an aperture size of about 5". 7 and an aperture correction factor of 1.73 
derived from the STinyTim PS F model. The absolute flux calibration uncertainty should 
be 4% (jEngelbracht et al.ll2006l ). Because of the contamination by the bright star a Ori in 
the IRAC and MIPS images, stars located within 30" of this O star were not included in the 
preliminary list. 



2.2. Optical photometry 



Optical magnitudes were obtained from the CIDA Equatorial V ariability survey which 
is being carried out using the QUEST I camera (IBaltay et al.l 120021 ) installed at the Jurgen 
Stock Telescope (a 1-m Schmidt telescope) at the Venezuelan National Observatory of Llano 
del Hato. The camera, an array of 4x4 CCDs, is designed to work in driftscan mode, which 
is a very efficient way to survey large areas of the sky. Each scan of the sky was reduced 



and ca librated with the standard QUEST software and the method described in I Vivas et al. 
(120041 ) for scans centered at 5=-l. We used data from multiple scans centered at 5 = — 3 



The scans have a width of 2.2° in 5. The range in right ascension shown in Figure 1 is covered 
completely by these optical observations. However, there is no data in a small region from 
-2.52° < 5 < -2.44°, corresponding to the gap between CCDs in the QUEST I camera. Most 
objects with an optical counterpart hav e more than 10 me asurements, so we can identify the 
variables using differential photometry (I Vivas et al.l 120041 ) . 



3. Membership selection 

We do not have spectroscopic confirmation for most of the probable members in the a 
Orionis cluster. We must use photometric criteria to make a robust selection of members, 
as described below. 
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3.1. Rejecting extragalactic objects 



Figure [3] shows the [3.6] versus [5.8]-[8.0] color-magnitude diagram for Sample 1 ( §2.1)) . 
The bulk of the sample has [5.8]-[8.0]~0 and [3.6]<14.5; below this limit the number of 
objects with photospheric colors ([5.8]-[8.0]~0) decreases dra stically, and the c ontamination 
from extragalactic sources is expected to be more than 50% (IFazio et al.l 12004 ). We identify 
extended sources (crosses) by computing the radial profile for each object using the radprof 
task in IRAF, and designate as extended sources those with FWHM values more than 3 
standard deviations from the median in all IRAC bands; typically this means FWHM>2". 
Most objects located below the [3.6] band limiting magnitude (solid line) in Figure [3] are 
identified in this way as extended sources and are likely to be extragalactic objects. We 
therefore initially restrict our sample to objects above the solid line in Figure [3j 

Some galaxies have very red colors because of strong PAH (polycyclic aromatic hydro- 
carbon) emissio n; with features at 7.7, 8.6, 11.2, 12. 7 , 16.4 , and 17.1 fim correlated with star 
formation rates ( Wu et al. 20051 ). Gutermuth et al. ( 2006 ) showed that the [4.5 - 5.8] versus 
[5.8 - 8.0] and [3.6 - 5.8] versus [4.5 - 8.0] color-color diagrams can be used to eliminate most 
of these objects. As shown in Figure HJ objects below and to the right of the solid lines are 
likely to be PAH-rich galaxies. We find that most of the PAH-rich galaxy candidates have 
extended radial profiles as inferred from point-source fitting. Visual inspection of the IRAC 
images also shows that many of these candidates are diffuse objects. In addition to the galax- 
ies, other objects with similar colors turn out to be faint companions of brighter stars, thus 
calling into question the accuracy of the photometry. Only one object (open square) with a 
stellar profile is located in the PAH-rich galaxies regio ns; this star (HD294268=SQ411 ) has 
been identified as a young star with spectral type F5 (jGregorio-Hetem fc Hetemll2002l ). We 
therefore eliminate the rest of the objects below the solid lines in Figure HI resulting in a 
sample of 1280 stellar sources (hereafter Sample 2). 

About 40% of the objects in Sample 2 do not have optical counterparts in the CIDA 
Equatorial Variability survey; most of these objects are fainter than the limiting magnitude 
(V~19.7) or brighter than the saturation limit (V~13.5) of this survey. Some of the objects 
in Sample 2 also are located on the gap i n declination of th e survey. We augmented the 
V optica l data set using pho tometry from (jSherry et al.l 120041 ) and for the brightest objects 
from the iKharchenkol (120011 ) catalog. 
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3.2. Additional photometric selection 



Virtually all the objects in Sample 2 (99%) have 2MASS counterparts fjCutri et al. 



20031 ). Figure [5] shows the distributions of J magnitudes for all 2MASS sources in the IRAC 
field (open histogram) and for the 2MASS sources included in Sample 2 (striped histogram). 
Comparison of these two histograms indicates that Sample 2 is essentially complete to ab out 
J=14.0. Assuming a reddening of E(B-V)=0.05 (Brown et al. 1994 ; Bejar et al. 19991 ) . a 
distance of D=440 p c (jBrown et al.lll998l ; iHernandez et al.l 120051 ) and the 3 Myr isochrones 



from lBaraffe et al.l (119981 ). the equivalent completeness limit in mass is ~ 0.15 M (spectral 
type ~ M5). For comparison, the distribution of J magnitudes for our MIPS detections 
(solid histogram) is shown in Figured Most of brightest objects (J<11; ~1.5 M Q ) have 24 
^m detections. 

Sample 2 includes 214 stars previously found to be members by other investigators 
(hereafter previously known members; see the Appendix |Aj). Of these objects, the member- 
ship of 104 stars were confirmed spectroscopically, while the other 110 stars were assigned 
membership based only on optical p hotometry. In ad d ition, 105 objects in our sample have 
been identified as X-ray sources by iFranciosini et al.l (120061 ). of which 49 objects have not 
been identified previously as spectroscopic or photometric members. 

Figure E] shows the V versus V-J color-magnitude diagram for Sample 2. Filled squares 
(spectroscopic members) and filled triangles (photometric members) represent previously 
known members; both distributions share the same region in Figure [6] suggesting that 
optical photometri c selection can be effecti ve in identifying low mass members (see also; 



Kenyon et al.ll2005l ; iBurningham et a 



the 3 Myr isochrone of iBaraffe et al. 



2005). The stars with 3 < V - J < 4 roughly follow 



(11998 



dashed curve), while the redder stars depart 
significantly from the theoretical isochrone; this may be due to incompleteness in the opac- 
i ty tables at low temperatures, resul ting in theoretical colors that are too blue for V-J > 4 



(ILvra et al.l 



2006 



Baraffe et al.lll998f) . Filled circles represent the higher mass stars selected 



by iKharchenko et al.l (120041 ) with high membership probabilities (>50%); the plus symbols 
represent bright stars with very low probabilities of membership. 

The mean of the distribution of previously known cluster members can be characterized 
roughly as straight line in the V versus V-J diagram, with a standard deviation cx(V-J) = 0.27 
mag. As a first approximation, we define the photometric region of probable members as the 
region above the dotted line in Figure El This line r epresents the lower 3 a deviation from the 
mean distribution line. The X-ray detected stars (iFranciosini et al.l 120061 ). labeled by open 
squares, are mostly located in the region of the previously known cluster members (most 
of these objects were identified previously as cluster members). Therefore, X-ray emission 
can be used as an additional criterion supporting membership in young stellar groups (e.g.; 
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Carkner et al.lll998l ; iPreibisch fc Feigelsonll2005l ; iFranciosini et al.ll2006l ) 



Figure [7J shows the J versus J-K color-magnitude diagram for Sample 2. Stars with 
J > 11.5 have similar J-K colors, in agreement with the trend expected for the 3 Myr 
isochrone from iBaraffe et al.l (119981 ). For these stars, the distribution of J-K color can be 
represented as a Gaussian function centered at J-K=0.94 with cr(J-K) = 0.06. We adopt the 
border between members and non- members to be the 3 a deviation line (dotted line). For 
stars with J<11.5, we define the region of probable members by using the X-ray sources in 
addition to the known members of the cluster. The crosses represent non-members defined 
using Figure [6] (below the dotted line). Some of these non-members, with J>13, appear 
in the region of probable members in Figure so some contamination would result from 
selecting stars using 2MASS colors alone. 

Our preliminary member s ample consists of stars located redward from dotted lines in 
Figures [6] and [3, X-ray sources (IFranciosini et al.ll2006l ). and previously known members (see 
£ |A]) located redward from the dotted line in Figure UJ Our preliminary uncertain member 
sample of the cluster includes those objects selected using only the criterion defined in Figure 
where the contamination of non-members could be important and additional confirmation 
is required. 

The most important source of contamination are field M dwarfs. Given the galactic 
latitude of the a Orionis cluster (b~17.3 deg), M gia nt stars are not expecte d to be a signif- 
icant source of contamination (<5%) in the cluster (IBejar et al.lll99S. 120041). We identify 3 
possible M giants stars using the J-H versus H-K color-color diagram (IBessell fc Brettlll9881 ). 
including them in the uncertain members sample. 



Jeffries et al.l ( 120061 ) showed that young stars located in the general region of the cluster 



consist of two populations with different mean ages and distances superimposed along the line 
of sight: one with an age of 3 Myr located at 439 pc (consistent with the a Orionis cluster) 
and the other wit h 10 Myr located a t 326 pc, which cou ld be associated with the Orion 
OBla association (jJeffries et al.l 120061 ; iBricefio et al.l 120061 ) . However, most of the objects 
belonging to the older stellar group are located north-ward from the main population of o 
Orionis cluster, so that the contamination by the older stars in our survey should be small, 
<10% (Jeffries 2006, private communication). In particu lar, 11 out 75 stars included in our 
preliminary member sample have radial velocities from iJeffries et al.l (120061 ) consistent for 
stars in the older population but inconsistent with those of the a Orionis cluster. These 
stars were included, with a note, in the uncertain member sample. Only two stars in our 
uncertain members sample have radial velocity expected of stars in the a Orionis cluster; we 
included these in our final membership list. 
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The final catalogs include 336 stars in the member sample and 133 stars in the uncertain 
member sample. In the member sample, 132 stars have optical photometry from the CIDA 
su rvey, of whic h 57 st ars (43%) were classified as variable stars using the method described 



in 



Vivas et al.l (120041); since y o ung stars are expected to exhibit photometric variability 



(e.g.; iCohen fc Schwartzlll976l ; iHerbst fc Shevchenkd Il999l ; iBriceno et al.l 120051 ) . this is an 
additional criterion supporting membership. 

Table [T] shows IRAC, MIPS and V photometry data for stars in the member sample 
according to the criteria discussed above. Column 1 shows the internal running identification 
number in our sample; columns 2 and 3 are the stellar coordinates; columns 4, 5, 6 and 7 give 
IRAC magnitudes in the bands [3.6], [4.5], [5.8] and [8.0], respectively; column 8 gives MIPS 
(24 jum) magnitudes; columns 9 and 10 show the V magnitudes and its sources; column 11 
indicates if the star has bee n identified as variable in the CIDA survey; column 12 indicates 
if the star is a x-ray source (IFranciosini et al.ll2006l ); column 13 shows the disk classification 
based on on the IRAC and MIPS analysis (see §1]); the last column gives the references 
for the previously known members. Table [2] shows IRAC and MIPS photometry for stars 
included in the uncertain member sample. The information shown in this table in columns 
from 1 to 8 is the same as in Table [TJ The last column shows the disk classification based 
on the IRAC and MIPS analysis discussed in the following section. 



4. Results 

4.1. Infrared color-color diagrams 

Figure [8] shows two IRAC color-color diagrams, [5.8]-[8.0] versus [3.6]-[4.5] (left panels) 
and [4.5]- [5.8] versus [3.6]- [4.5] (right panels), for the members (upper panels) and uncertain 
members (lower panels). Different symbols show the disk classification based on the IRAC 
SED slope ( §4.2p . Particularly, triangles, squares and diamonds represent stars with different 
levels of infrared emission, while circles are non-excess stars. Symbols surrounded by open 
squares show objects with small or no excesses in the IRAC bands but large excesses at 24 
/im (see below). The dashed box in the [3.6] - [4.5] versus [5.8] - [8.0] color-color diagram 
indicates the th e region predicted for stars with optically-thick disks ( §4.2p by the models of 



D'Alessio et all fl2005ri ) 



The upper-left panel of Figure [8] shows that most of stars with significant IRAC excess 
(squares; hereafter IRAC class II stars) have colors in agreement with model colors for 
optically-thick disks. The stars without IRAC excesses (hereafter IRAC class III stars) are 
located around the photospheric region ([5.8]-[8.0]~0, [3.5]-[4.5]~0). A few systems are found 
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in the region between the IRAC class III and class II stars; these objects with weak IRAC 
excess (triangles) are discussed in more detail in §4.21 Two of these objects (the stars S0638 
and S0299) have [3.6]- [4.5] <0, while exhibiting significant emission at longer wavelengths. 
Three other objects (diamonds) have very large excesses which may suggest the presence of a 
protostellar envelope (he reafter IRAC class I cand idates); S0457 (IRAS05358-0238), S0927 



an emission line object IWeaver fc Babcockl l2004f ). and S0411 (HD294268). The last two 



objects have [3.6]-[4.5] colors expected for class II systems, so these objects are more likely 
to be objects with optically-thick disks and high 8 /xm fluxes due to strong silicate emission 
or rising SEDs at wavelength longer than 5.8 jum ( §4.3p . 

In the upper-right panel of Figured], the [4.5]-[5.8] colors of IRAC class II stars shows 
less scatter in comparison with the [5.8]-[8.0] colors. There may be some PAH background 
contamination present in the [8.0] band; however, this also could indicate that intrinsic 
differences in disk fla ring or other properties that may be more easily detected at 8/iin 



(ID'Alessio et al.l 120061 ). For example, the IRAC class II stars located redward from the 
dashed box (S0462, S0774 and S01266) and the star S0411 (with [5.8]-[8.0]>1.5) in the 
upper-left panel, have a flat or rising SEDs at wavelengths longer than 5.8 /j,m (see §4.31) . 

In the lower panels, nine IRAC class II stars are located in the region expected for stars 
with optically thick disks, and thus are likely to be cluster members. One IRAC class II star 
(S01356), located below the model region, has an excess at [5.8]-[8.0], but not at [3.6]-4.5] 
or [4.5]-[5.8]. This star was not covered in the MIPS field, so it is not clear if the excess 
observed at 8 /im originates from PAH background contamination or from emission from an 
outer disk. The object S01293 (represented by the triangle) has a larger excess at [4.5]- [5.8] 
than at [5.8]-[8.0]; it is faint (J=15.64) and so could have a longer wavelength excess below 
the MIPS detection limit. 

Five uncertain members with very large infrared excesses (IRAC class I candidates), are 
plotted in the IRAC color-color diagrams: SO950, S0336, S0361, S0916 and SQ6 68. Because 



AGN are located in the same region as the class I objects in these diagrams (jStern et al. 



20051 ) and are too faint to easily detect spatial extension, additional data are needed to 



confirm these objects as protostellar members. 

Figure [9] shows K-[24] versus J-H color-color diagrams for members (upper panel) and 
uncertain members (lower panel). We sort th e objects into three regi ons, depending upon the 



magnitude of their K-[24] excesses (see also iHernandez et al.ll2006l ). One region consists of 



K-[24] < 0.5, corresponding to systems with undetectable or very uncertain long-wavelength 
excesses, given photometric errors. A second region encompasses objects with K-[24] >3.3, 
corresponding to optically-thick disk emission and/or protostellar envelope emission. Finally, 
the intermediate-excess region 0.5< K-[24] < 3.3 is consistent with colors typical of optically- 
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thin debris disks [the famous f3 Pic debris disk system, whic h is one of the debris disk system 



with largest 24 /im excess known, has K-[24]~3.3; (e.g., iGorlova et al.l l2004j ; iRieke et al. 



20051 )]. 



There are seven members with very large [24] excesses but low IRAC excesses (symbols 
surrounded by open squares); SO908, S0818, S01267, S0897, S0299, S0587 and S01268. 
The last object has no detectable IRAC excess. These stars are possible "transition disk" 
systems, in which the inner disk is c ompletely or nearl y-completely cleared of small dust 
but the outer disk is optically thick (ICalvet et al.ll2002l . see §4.11) . In the lower panel, the 
uncertain member SO120 is also a transition disk candidate. 

There are two early type members which are located in the intermediate-excess region: 
S0913 (J-H=-0.04) has no IRAC excess and S0956 (J-H=0.10) has a weak IRAC excess 
( §4.2p . These stars are "debris disk candidates". The star S0981 (J-H~0.45) shows no 
IRAC excess but has a modest 24 /zm excess, and thus is a likely debris disk system around 
a K early type star. However, a more detailed study of the debris disk candidates in this 
young stellar cluster is necessary to clarify the actual nature of their 24 /xm excess es, since 
contribution from primordial disk material could be present ([Hernandez et al.ll2006l ). There 
are 5 low mass star members (J-H>0.5) with weak IRAC excess, in which the excesses 
observed at the IRAC bands and at 24 fim are consistent with a overall decrease of infrared 
emission from the inner and outer disks. 



4.2. IRAC SED slope 



Lada et al.l (120061 ) showed the utility of characterizing excess infrared emission using the 
IRAC SED slope a determined from the [3.6]-[8.0] color, where a=dlog[\F\]/dlog[\]). In 
this system, IRAC class I candidates, defined as systems with a > 0, are usually protostars; 
disk-bearing systems, having —2.56 < a < 0, include IRAC class II stars and objects with 
weak IRAC excess (the lower limit depends upon spectral type and photometric errors, see 
below); and IRAC class III stars have essentially no inner disk emission, with a < —2.56 
(again, depending upon spectral type and errors). 

Figure [10] shows a versus the [8.0] magnitude for members (filled symbols) and uncertain 
members (open symbols) of the cluster. Dashed lines represent the regions of differing levels 
of infrared excess as discussed above. We subdivide the disk-bearing stars into two categories: 
IRAC class II sta rs (squares, a>-1.8); and systems with less IRAC infrared excess, which 
Lada et al.l (120061 ) termed "anemic" disks and which we call here "evolved disks" (-2.56 



< a < -I. 
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Because the photometric uncertainties in a increase with [8.0], a single straight line 
in Figure [10] to sepa rate IRAC c l ass H I stars (in the IRAC wavelength range) and evolved 
disks is not useful. lLada et al.l (120061 ) avoided this problem by applying their criteria to 
separate class III and evolved disks only to the brighter stars with spectral types M4.5 or 
earlier (the number of evolved disks would be overestimated in the lowest mass bins due to 
large photometric errors; C. Lada, private communication). We distinguish between IRAC 
class III stars (circles) and evolved disk systems (triangles) using the mean error of a (a) 
calculated from the photometric errors for stars within 0.5 magnitude bins. IRAC class III 
stars are those with excesses, if any, smaller than the 3 a deviation (solid lines) from the 
median a. 

We wish to examine the mass dependence of disk emission, but we do not have spectral 
types for most of these stars. As a proxy for mass, we use the J magnitude, which is relatively 
unaffected by disk excess emission, and convert to approximate mass ranges using the 3 Myr 
isochrone from lBaraffe et al.l (119981 ) for low mass (<1 M©) stars and lSiess et al.l (120001 ) high 
mass (>1 M0) stars, assuming a distance of 440 pc and E(B-V)=0.05. Figure [TT1 shows a 
versus J, with the regions expected for objects in different ranges of mass: HAeBe range (>2 
M©), Intermediate Mass TTS (IMTTS) range (1.0-2.0 M©), TTS range (0.1-1.0 M©) and 
Brown dwarfs range (<0.1 Mq). Most of the objects with optically thick and evolved disks 
are located in the mass range of T Tauri stars. 

Using the results from Figure [TT], we derive the fractions of stars with optically-thick 
disks (filled circles + dotted line) and all disk-bearing stars (open circles + dashed line) 
versus J magnitude, with corresponding mass ranges indicated (see also Table [5]). The lowest 
fraction of disks is observed in the highest mass stars (HAeBe range), i n agreement with 



the fr action calculated for the overall population of OBlb sub-association ([Hernandez et al. 



20051 ). There is a marginal evidence that the disk fraction de clines towar d lowe r masses 
(into the brown dwarf range), in agreement with results from lLada et al.l (120061 ) for the 
young (2-3 Myr) stellar group IC 348. Nevertheless, including the error bars, the fraction is 
also consistent with no mass dependence toward l ower masses. The fra ction of thick disks 
in the TTS mass range agrees with the result from lOliveira et al.l (120061 . 33±6 %; ) for low 
mass stars. 

Only 15±5% of disk-bearing stars in the TTS mass range are disk evol ved objects, ap- 
proxi mately in agreement with results from the 4 Myr-old cluster Tr 37 (10%; ISicilia-Aguilar et al 
20061 ). although with better statistics and higher photometric accuracy. If we assume a 
steady-state evolution of disks, and also assume that the evolved disks are going through 
a phase of clearing of small dust particles in the inner disk, then our observed frequency 
suggests a timescale of inner disk clearing ~ 0.15 x 3 Myr ~ 0.45 Myr. 
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4.3. Disk diversity: SEDs 



Figure QJ shows SEDs (data from IRAC, MIPS, 2MASS and optical B, V, R and I 
magnitudes when are available) for selected members. The first 5 rows show the SEDs of 
stars in the mass range of TTS (0.1-1.0 M©), while the last row shows the SEDs of stars in 
the mass range of HAeBe stars (>2 Mo). SEDs (solid lines) are normalized to the J band and 
are sorted by different levels of IR excesses indicating differences in the disks. Each panel is 
labeled with the internal running identification number and the classification from the IRAC 
SED slope analysis ( §4.2p . The dotted line represents the median SED (Table [3]) of optically 
thick disk stars in the mass range of TTS defined in Figure [TT1 (hereafter, optically thick disk 
median SED); error bars denote the quartiles of the distri bution. This typical d isk emission 
in the a Orionis cluster is lower than the median in Taurus (ID'Alessio et al.lll999l . dot-dashed 
line in the first upper-left panel), which has an age of ~l-2 Myr, indicating differences due 
to evolutionary effects such as grain growth and/or settling to the disk midplane (see als o 
Sicilia-Aguilar et al.ll2006l ; lLada et al.ll2006l ). or lower accretion-rate (ID'Alessio et al.ll2006l ). 
The dashed lines represent the median photospheric SED (Table H|) in the T Tauri mass 
range. 



The first row of Figure [13] shows the SEDs of evolved disk objects (SO905, SO 1057, 
S0759 and S0587) which exhibit modest infrared disk emission. The lower level of infrared 
excess emission could be explained by decreasing the height of the irradiation surface in the 
disks due to a higher degr ee of settling, which d iminishes the degree of flaring, thus produc- 
ing flatter disk structures Jp'Alessio et al.ll2006h . The second row of SEDs (SO120, S01268, 
S0299 and S0897) shows some transition disk candidates, which exhibit little or no excess in 
IRAC bands but 24 /zm emission similar to optically thick disks. This indicates a removal of 
small dust particles in the inner disks but not depletion or removal in the outer disk. In the 
next rows, the stars S0818, S01267, SO908, S01156, SO540 and S01266 show somewhat 
weaker IRAC excesses. The stars S073, S0927, S01154, SOI 153, S0457, and S0668 show 
larger infrared excess es than the optically t hick median, which might in dicate more highly 
flared disk structure (ID'Alessio et al.ll2006l ; iDullemond fc Dominikl 12005), edge-on disks, or 
a con tribution from an envelope. S0927 has strong emission in the Ha ( I Weaver fc Babcock 
20041 ). S01153 (V510 Ori; HH 444) and SO 1154 (HH 445) are Herbig-Haro objects pro- 
ducing collimated outflows ([Andrews et al.l 120041 ) . Finally, S0457 has been identified by 
Oliveira &: van Loonl (120041 ) as a class I object. 



In the last row of Figure Q31 SEDs of earlier-type membe rs are compared with photo- 
spheres of stars with similar spectr al types (long dashed line iKenyon fc Hartmannlll995l ). 
The star S0724 (B2; lGuetterlll98l[ ) shows small excess at 8 /xm but has no excess at 24 /xm, 
this star was classified in §4.21 as an evolved disk object possibly due to a background PAH 
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cont a nimation at 8 /mi. The stars S0913 (B9.5; iGuetterl Il98lh and S0956 (A8; iGuetter 
198ll ) show small or no excesses in the IRAC bands but have significant excesses at 24 /im, 
consistent with being debris disk stars. Finally, the star S0411, classified as class I candi- 
date in §4.21 shows small excesses at wavelengths shorter than 5. 8 /mi and a rising SEP at 
longe r wavelengths. This emission line star, with spectral type F5 (jGregorio-Hetem fc Hetem 
20021 ) . could be an accreting transition disk object in the HAeBe mass range. This object 
alternatively could be considered t o be an IMTTS, as F5 i s the usual transition spectral type 



between HAeBe and IMTTS (see IHernandez et al.ll2004l ). 



Most of the disk classifications using the IRAC SED slope ( §4.2j) agree with the results 
from analyzing SEDs. However, a few objects were reclassified taking into account the overall 
SEDs; these objects are labeled in Tables [JJ and [2J 



4.4. Disk evolution 



Figure [TH shows the fraction of stars in the TTS mass range with near - infrared disk 



emiss ion in different stellar groups, as a function of age (IHernandez et al.ll2005l; lHaisch et al. 



2001). Recent Spitzer results 
20041 ). Ta urus (IHartmann et a 



NGC7160 flSicilia-Aguilar et al. 



or the young stellar clus ters NGC7129 (iGutermuth et al. 



2005ch . 7] C hameleontis dMeeeath et ali l2005ah. Tr37 and 



2006|), IC348 JLada et al-kood ). and Upper Scorpius ([Carpenter et aJj 



20061 ) have been included in this plot. The disk fractions shown in Figure [141 include stars with 
optically thick and evolved disks (disk-bearing stars) detected mostly using IRAC photome- 
try. The a Orionis cluster (filled square) follows (with in the scatter of the distribution) the 



disk fraction trend outlined by in other stellar groups (lHaisch et al.ll200ll ; iHillenbrand et al. 
2006h . 



Aside from the decrease in the disk fraction with age, the amount of infrared disk 
emission also decreases wi th age. Figure [T5l shows the dist ribution of a for l ow mass stars in 
the stellar groups Taurus ([rlartmann et al. 2005c ). IC348 ( Lada et al. 2006 ). o Orionis (this 



work) and Tr 37 (jSicilia-Aguilar et al 



20061 ). As in Figure [TBI the youngest group, Taurus, 



has a larger population of stars with robust optically thick disks [a > -1). In contrast, in the 
more evolved stellar groups most disk-bearing stars have an a <-l, indicating a reduction in 
the height of the disk photosphere (less disk flaring) in the inner regions. 

The Taurus sample suggests a rel atively clear break betw een stars without IRAC ex- 
cesses and optically-thick disk systems (IHartmann et al.ll2005a ). although the sample size is 
relatively small. The a Orionis sample also suggests a small fraction of objects with inter- 
mediate or evolved disk excesses, as discussed previously. The evolved disk fraction is less 
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certain in Tr 37, due to the larger photometric errors for that more distant cluster. The IC 
348 sample suggests a larger fraction of evolved disks than seen in a Ori; the reality of this 
possible difference requires further study. 

In contrast to IC 348, the a Orionis and Tr 37 clusters include high mass stars located 
in their center (09 in a Orionis cluster and 07 in Tr37); however, the medians of a are 
similar in these stellar groups (Figure [T5l) . This suggests that external photoevaporation in 
o Orionis and Tr 37 may be not very important for disk dissipation. 



4.5. Spatial distribution of the a Orionis cluster 



Figure [161 shows the spatial distribution of members with and without infrared excesses. 
As shown in the upper right panel, the space density has a peak with a FWHM of 0.25° 
(1.9 pc). The upper-left panel shows the ratio of systems with disks to the total members 
in radial distance bins defined by the semi-circles East-West versus the radial distance from 
the center of the cluster. Solid circles (+dotted line) represent the fraction of thick disks, 
while open circles (+dashed line) represent the total fraction of disk-bearing stars, defined 
in §4.31 The error bars represents the statistical N~ l l 2 error in our disk fractions, where N 
is the number of stars with disks. Similarly, the lower-right panel shows the fraction of disks 
in radial distance bi ns defined by the semi-circles North-South. In contrast to the result of 
Oliveira et al.l (120061 ). using K and L band observations, there is evidence for a higher disk 
fraction near the cluster center. However, we cannot rule out a small, relatively uniform 
contamination of non-members which could also produce this effect. 

There is a slight increase in the disk fract ion in the north-ea s t direction; in the directio n 
toward NGC2024 (age < 1 Myr, D~440 pc; iLevine et allbood : lEisner fc Carpenter! l2003h . 
located at ~1.0 deg from the a Orionis center; this fraction could be affected by disk-bearing 
members of this younger stellar group. Since NGC2024 is centrally concentred, the fraction 
could be also affected by stars extends toward the Horsehead nebula (located at ~0.6 deg 
from t he a Orionis center), which also shows several class II objects in separate Spitzer 
maps (IMegeath et al.ll2005bl ). Spectroscopic studies are necessary to test whether the disk 
fractions in the outer regions are affected by non-members. 



4.6. Disk and X-ray emission 



A sub-sample of 81 stars have X-ray luminosities from lFranciosini et al.l (120061 ). including 
8 evolved disks, 21 optically thick disks and 52 class III stars. Figure [TTl shows the cumulative 
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distribution of class III stars and stars with optically thick disks as a function of X-ray 
luminosity. The significance level of 64% in a Kolmogorov-Smirnov (KS) test indicates that 
class III stars have similar X-ray luminosities as stars with disks. This suggests that the 
presence of optically thick disks have small or no influence on X-ra y emission, such a s migh t 
be caused by an accretion shock, in agreement with the results from lFranciosini et al.l (120061 ). 



4.7. Contamination expected in the uncertain members sample 

Stars selected as uncertain members are relatively faint (J>13.0). From 133 stars in 
this sample, only 11 objects have IRAC SED slopes consistent with the presence of disks 
(10 thick disk + 1 evolved disks). So, the fraction of disk-bearing stars for the uncertain 
members (8.3±2.5 %) is very low compared with the fraction of disk-bearing stars expected 
for members with J>13 in the a Orionis cluster (31.3±4.8). Assuming that the source of 
non-members are class III stars and the disk-bearing stars are actual members of the cluster, 
we can expect that at least 84 class III stars (~70 % of the uncertain member sample) do 
not belong to the a Orionis cluster and are likely to be field M dwarfs or other sources of 
contaminations. Additional spectroscopic studies are essential to determine membership for 
this sample. 



5. Conclusions 

We have used the IRAC and MIPS instruments on board the Spitzer Space Telescope 
to study the frequencies and properties of disks around 336 members of the 3 Myr old a 
Orionis cluster, with results for another 133 objects of uncertain membership. We confirm 
and strongly improve upon previous results showing that disk fractions are much lower 
around A-type and intermediate-mass T Tauri stars, implying that disk evolution is much 
more rapid for stars with masses > 1M Q . In particular, we find that only a few percent 
of all A-F stars have optically thick disks, with perhaps another ~ 15 % having detectable 
optically-thin, probable debris disks. This fraction of probable debris disk is lower that those 
found in the older Orion OBlb (38%) and OBla (46%) sub-associations, suggesting that cr 
Orionis cluster could be in a phase between the clearing of the primordial disk and the 
formation of l arge icy objects, which produce the second-generation dust observed in debris 



disk systems ( IHernandez et al.ll2006l ). The disk frequency ~ 30 — 35% we find for T Tauri 
stars is reasonably consistent with studies of other populations, perhaps a bit lower. We also 
find no evidence that the disk fraction decreases significantly among the highest-mass brown 
dwarfs. Overall, disk frequencies follows the overall trend of decreasing disk frequency with 
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increasing age observed in other young stellar populations. 

About 15% of the T Tauri disk systems exhibit weak IRAC excesses inconsistent with 
simple optically-thick disk models. If these systems are in the process of clearing small dust 
particles from their inner disks, the observed fraction suggests that the timescale for this 
evolution is < 0.5 Myr. We also find that the infrared excesses among optically-thick disks 
in a Ori are systematically lower than the excesses observed from younger Taurus molecular 
cloud stars. This result confirms and strongly extends previous results indicating an overall 
decrease in the level of emission of optically-thick disks with age, which may indicate an 
overall settling of dust particles to the disk midplanes, or grain growth, or both. The 
populous nature and relatively near distance of the a Orionis cluster allows us to improve 
the statistics of these correlations, providing an important step toward a full characterization 
of protoplanetary disk evolution in clustered environments, in which most stars form. Finally, 
we do not find correlation between X ray luminosity and the disk infrared excess, suggesting 
that the presence of optically thick disks have small or not influence on the X ray emission. 

Thanks to R.D. Jeffries, D. Stern and C. Lada for insightful communications. This pub- 
lication makes use of data products from the CIDA Equatorial Variability Survey, obtained 
with he J. Stock telescope at the Venezuela National Astronomical Observatory, which is 
operated by CIDA for the Ministerio de Ciencia y Tecnologia of Venezuela, and Two Micron 
All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared 
Processing and Analysis Center/ California Institute of Technology. This work is based on 
observations made with the Spitzer Space Telescope (GO-1 0037 and GO-1 0058), which is 
operated by the Jet Propulsion Laboratory, California Institute of Technology under a con- 
tract with NASA. Support for this work was provided by University of Michigan and NASA 
Grant NAG5-13210. 



A. Previous membership studies 



The cr Orionis cluster has been the subject of several photometric and spectroscopic 
studies to identify members of its population.. A nalyzing the photometric and kinematic 
properties of the bright population in the cluster, iKharchenko et al.l (120041 ) determined an 
angular size of 0.40 degrees and identified 7 stars with photometric and kinematic member- 
ship probabilities larger than 50%. Since the cluster is reasonably close, relatively young 
and has low reddening, many of the previous studies have focused on identifying low mass 
stars and substellar objects. Adopti ng an age in t erval of 2-7 Myr and the Hipparcos dis- 
tance of the star a Ori (352i^ 6 pc), Ifjeiar et all Jl999l . I2OO1L E3) presented several lists 
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of low mass stars, brown dwarfs and free-floating planets belonging to the a Orionis cluster, 
selected from optical and NIR color-magnitude diagrams. Adopting an age o f 2.5 Myr and 
the distance calculated for OBlb subassociation (440 pc). ISherry et all (120041 ) reported pho- 
tometric membership for a sample of low mass stars (193 objects) based on the V versus V-I 
color-magnitude diagram. ISherry et al.l (120041 ) found that the cluster has ~160 members in 
the mass range from 0.2 M to 1.0 M Q . Adopting an age of 3 Myr and a distance of 350 pc, 
Scholz fc Eisloffell (120041 ) reported photometric membership for 135 candidates of very low 
mass stars (<0.4 M ) using the color magnitude diagrams, I versus I-J, and photometric 
variability. 



Zapatero Osorio et al.l (120021 ) presented intermediate and low resolution spectra includ- 



ing the Ha and Li I A6 708A lines and confirmed 25 low m ass stars and 2 brown dwarfs as 
members of the cluster. iBarrado y Navascues et al.l (120031 ) confirmed 51 low mass stars and 
brown dwarfs as members using l ow resolution s pectr oscopic analysis of forbidden lines, as 
well as Li I A6708A and Ha lines. iKenyon et al.l (120051 ) presented an intermediate-resolution 
spectroscopic study of 76 photometrically selected stars ranging in mass from 0.04 M Q to 
0.3 M ; 57 stars were confirmed as members on the basis that they exhibit Li I A6708A 
in absorption, a weak Na I doublet (8183,8195 A) indicating low surface gravity consistent 
with PMS status, a nd ra dial velocities expected for the kinematic properties of the cluster. 
Burningham et al.l (120051 ) confirmed as members 38 low-m ass s tars using high-res o lution 
spectroscopic analysis of Na I doublet. IKenyon et al.l (120051 ) and iBurningham et al.l (120051 ) 
showed that optical photometric selection alone is very effective to identify low mass mem- 
bers belonging to the cluster, with small contamination by non-members for stars with I<17 
magnitudes. 



Recently, lOliveira et al.l (120061 ) compiled 59 spectroscopic confirmed members of the 



cluste r and studied the circumstellar disks of the sample using NIR data. iFranciosini et al. 



(120061 ) detect 175 X-ray sources down to the substellar limit. Half of this sample could be 
new members of the cluster based on the X-ray emission as indicator of youth. 
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Table 1. Members of the a Orionis Cluster 



so 


RA(2000) 


DEC(2000) 


[3.6] 


[4.5] 


[5.8] 


[8.0] 


[24.0] 


V 


Rcf 


Var 


X-ray 


Disk 


Known 


ID 


dcg 


dcg 


mag 


mag 


mag 


mag 


mag 


mag 






source 


type 


member 



9 


84.32781 


-2 


67275 


10.99 ± 0.03 


10.97 ± 0.03 


10.80 ± 0.03 


10.92 ± 0.03 








15.40 


± 0.05 


q 







III 9 


27 


84.34599 


-2 


54623 


13.00 ± 0.03 


12.94 ± 0.03 


12.86 ± 0.04 


12.99 ± 0.05 








18.42 


± 0.02 


9 






III 


41 


84.35401 


-2 


86730 


10.30 ± 0.03 


10.31 ± 0.03 


10.28 ± 0.03 


10.25 ± 0.03 


















III 


52 


84.36209 


-2 


36497 


10.54 ± 0.03 


10.56 ± 0.03 


10.50 ± 0.03 


10.49 ± 0.03 








13.77 


± 0.08 


q 







III 


59 


84.36683 


-2 


60175 


12.65 ± 0.03 


12.62 ± 0.03 


12.51 ± 0.04 


12.52 ± 0.04 








17.64 


± 0.05 


q 







III 


60 


84.36792 


-2 


40508 


12.87 ± 0.03 


12.82 ± 0.03 


12.82 ± 0.04 


12.80 ± 0.05 








18.08 


± 0.11 


q 







III 9 


73 


84.37888 


-2 


39521 


10.90 ± 0.03 


10.44 ± 0.03 


10.04 ± 0.03 


9.11 ± 0.03 


5 


(14 


± 0.03 


16.28 


± 0.10 


q 


1 




II 9 


74 


84.37941 


-2 


52922 


9.46 ± 0.03 


9.39 ± 0.03 


9.30 ± 0.03 


9.32 ± 0.03 


9 


19 


± 0.04 












III 


77 


84.38132 


-2 


40748 


10.99 ± 0.03 


10.98 ± 0.03 


10.93 ± 0.03 


10.93 ± 0.03 








15.83 


± 0.11 


q 


1 




III 9 


82 


84.38269 


-2 


75514 


9.24 ± 0.03 


9.21 ± 0.03 


9.16 ± 0.03 


9.20 ± 0.03 


8 


95 


± 0.04 


9.42 ± 0.02 








III 


92 


84.38811 


-2 


89332 


10.59 ± 0.03 


10.74 ± 0.03 


10.64 ± 0.03 


10.55 ± 0.03 








14.99 


± 0.04 


q 







III 9 


115 


84.40181 


-2 


24377 


10.74 ± 0.03 


10.72 ± 0.03 


10. 07 ± 0.03 


10.70 ± 0.03 








15.49 


± 0.10 


q 







III 9 


116 


84.40192 


-2 


69906 


14.19 ± 0.03 


14.02 ± 0.03 


13.98 ± 0.06 


13.97 ± 0.09 


















III 1,2,3,6 


117 


84.40270 


-2 


56669 


11.82 ± 0.03 


11.83 ± 0.03 


11.81 ± 0.03 


11.76 ± 0.04 








17.12 


± 0.05 


q 







III 9 


123 


84.40760 


-2 


76226 


11.51 ± 0.03 


11.47 ± 0.03 


11.40 ± 0.03 


11.44 ± 0.03 








16.37 


± 0.05 


q 


1 




III 9 


124 


84.40766 


-2 


83598 


9.49 ± 0.03 


9.54 ± 0.03 


9.45 ± 0.03 


9.44 ± 0.03 


9 


(15 


± 0.05 












III 


147 


84.42406 


-2 


48557 


9.51 ± 0.03 


9.50 ± 0.03 


9.47 ± 0.03 


9.53 ± 0.03 


9 


59 ± 0.04 


10.90 


± 0.07 


7 






III 


164 


84.43706 


-2 


49920 


9.91 ± 0.03 


9.90 ± 0.03 


9.91 ± 0.03 


9.89 ± 0.03 


9 


80 ± 0.05 


10.95 


± 0.09 


7 






III 


165 


84.43856 


-2 


48111 


12.65 ± 0.03 


12.58 ± 0.03 


12.55 ± 0.04 


12.49 ± 0.04 








18.76 


± 0.13 


9 






III 9 


168 


84.43889 


-2 


73678 


8.96 ± 0.03 


9.02 ± 0.03 


8.95 ± 0.03 


8.95 ± 0.03 


8 


94 ± 0.04 


10.66 


± 0.07 


7 






III 


194 


84.45671 


-2 


60506 


10.12 ± 0.03 


10.17 ± 0.03 


10.06 ± 0.03 


10.10 ± 0.03 


10.03 ± 0.06 












III 


2(19 


84.46288 


-2 


43542 


13.64 ± 0.03 


13.53 ± 0.03 


13.68 ± 0.05 


13.38 ± 0.07 


















III 2,3 


214 


84.46497 


-2 


59045 


10.81 ± 0.03 


10.80 ± 0.03 


10.79 ± 0.03 


10.72 ± 0.03 








15.64 


± 0.01 


9 




1 


III 8,9 


219 


84.46685 


-2 


60126 


13.92 ± 0.03 


13.81 ± 0.03 


13.73 ± 0.05 


13.70 ± 0.07 



















III 6 


220 


84.46747 


-2 


56052 


11.85 ± 0.03 


11.85 ± 0.03 


11.78 ± 0.03 


11.77 ± 0.04 








16.89 


± 0.03 


9 







III 9,12 


22!) 


84.47091 


-2 


55951 


9.45 ± 0.03 


9.46 ± 0.03 


9.43 ± 0.03 


9.42 ± 0.03 


9 


40 


± 0.04 










1 




240 


84.47515 


-2 


74462 


11.92 ± 0.03 


11.88 ± 0.03 


11.86 ± 0.03 


11.82 ± 0.03 








17.02 


± 0.05 


q 







III 8,9 


243 


84.47665 


-2 


65823 


8.51 ± 0.03 


8.68 ± 0.03 


8.55 ± 0.03 


8.51 ± 0.03 


8 


47 


± 0.03 


10.91 


± 0.07 


7 




1 


III 7 


244 


84.47683 


-2 


72711 


8.99 ± 0.03 


9.06 ± 0.03 


8.98 ± 0.03 


8.96 ± 0.03 


8 


98 


± 0.04 














247 


84.47852 


-2 


68587 


12.04 ± 0.03 


11.80 ± 0.03 


11.45 ± 0.03 


10.82 ± 0.03 


8 


41 


± 0.03 


18.55 


± 0.06 


q 








II 5,9 


251 


84.47960 


-2 


46006 


12.45 ± 0.03 


12.36 ± 0.03 


12.30 ± 0.04 


12.46 ± 0.04 








18.49 


± 0.10 


9 






III 9 


254 


84.48156 


-2 


55146 


13.66 ± 0.03 


13.35 ± 0.03 


12.94 ± 0.04 


12.32 ± 0.04 


9 


70 


± 0.05 













II 2,3,12 


256 


84.48218 


-2 


40933 


13.79 ± 0.03 


13.74 ± 0.03 


13.56 ± 0.06 


13.77 ± 0.10 


















III 2,3 


2(10 


84.48530 


-2 


73263 


11.91 ± 0.03 


11.96 ± 0.03 


11.85 ± 0.03 


11.86 ± 0.03 








16.84 


± 0.05 


q 







III 9 


271 


84.48935 


-2 


64562 


12.92 ± 0.03 


12.76 ± 0.03 


12.58 ± 0.04 


12.03 ± 0.04 


9 


40 


± 0.04 













II 1,2,3,6,8,12 


275 


84.49076 


-2 


91668 


11.62 ± 0.03 


11.62 ± 0.03 


11.60 ± 0.03 


11.55 ± 0.03 








14.86 


± 0.04 


q 


1 






283 


84.49328 


-2 


69059 


12.06 ± 0.03 


11.96 ± 0.03 


11.90 ± 0.03 


11.96 ± 0.04 








18.92 


± 0.07 


q 








III 6,9,12 


293 


84.49872 


-2 


85090 


9.61 ± 0.03 


9.70 ± 0.03 


9.62 ± 0.03 


9.59 ± 0.03 


9 


(17 


± 0.04 














297 


84.50228 


-2 


75267 


11.56 ± 0.03 


11.47 ± 0.03 


11.47 ± 0.03 


11.43 ± 0.03 








17.75 


± 0.06 


q 





1 


III 9,12 


299 


84.50394 


-2 


43548 


11.68 ± 0.03 


11.76 ± 0.03 


11.66 ± 0.03 


11.33 ± 0.03 


6 


7(1 


± 0.03 


16.70 


± 0.11 


q 


1 




TD 2 9 


300 


84.50434 


-2 


76042 


10.12 ± 0.03 


9.70 ± 0.03 


9.29 ± 0.03 


8.64 ± 0.03 


5 


87 


± 0.03 


17.59 


± 0.07 


q 


1 




II 9 



Note. — Table [Tl is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding its form and content. 

References in col um ns 10 and 14: q CIDA su rve y; 1 iBar rado v N avaaeueS et alj 12003 1); 2l5eiar et al.l 11999ft: SlBemr et al.l 120011): 4 IB6iar et all 12004(1; 5lBurn .ingh.am et all 120051) ; 6 
iKenvon et al.l 120051) ; 7 Kharchcnko ct al ] 120041) ; 8 lQliveira et al.l 120061 )7 120041) ; lO lScholz & Eisloffell 120041) ; ll fZapatero Osorio et al.l 120021) ; l^ Tjeffrie^^t^il' 1 H2006T ) 

Column 11: 1 variable; non-variable 



Column 12: 1 X-ray emission; not X-ray emission 

Column 13: III— non-excess, II — thick disk, I — class I candidate, EV— evolved disk, TD— transition disk candidate, DD — Debris disk candidate 

1 Classified as III in g3] 

2 Classified as EV in i|4"2l 
Classified as II in JOl 

4 Classified as I in JOl 



OX 
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Table 2. Uncertain members of the a Orionis Cluster 



so 


RA(2000) 


DEC (2000) 


[3.6] 






[4.5] 






[5.8] 






[8 


.0] 






[24.0] 


Disk 


ID 


deg 


deg 


mag 






mag 






mag 






mag 






mag 


type 


33 


84.34808 


-2.74109 


14.49 ± 





03 


14.46 ± 





04 


14.42 


± 





07 


14.44 


± 





15 




111 


17 


84.35836 


-2.34652 


13.53 ± 





03 


13.55 ± 





03 


13.47 


± 





05 


13.80 


± 





11 




111 


57 


84.36604 


-2.51190 


13.45 ± 





03 


13.45 ± 





03 


13.33 


± 





04 


13.52 


± 





08 




ill 


71 


84.37812 


-2.40533 


14.45 ± 





03 


14.48 ± 





04 


14.51 


± 





11 


14.21 


± 





13 




111 


81 


84.38234 


-2.40128 


14.36 ± 





03 


14.27 ± 





04 


14.19 


± 





07 


14.13 


± 





11 




111 


87 


84.38534 


-2.75422 


13.67 ± 





03 


13.70 ± 





04 


13.86 


± 





08 


13.96 


± 





15 




111 


107 


84.39636 


-2.44933 


12.10 ± 





03 


12.01 ± 





03 


11.91 


± 





03 


11.94 


± 





03 




III 


120 


84.40674 


-2.42888 


13.48 ± 





03 


13.49 ± 





03 


13.35 


± 





05 


13.30 


± 





09 


9.95 ± 0.05 


TD 1 


219 


84.46685 


-2.60126 


13.92 ± 





03 


13.81 ± 





03 


13.73 


± 





05 


13.70 


± 





07 




III 6 


127 


84.40806 


-2.80576 


12.78 ± 





03 


12.82 ± 





03 


12.80 


± 





04 


12.73 


± 





04 




III 


134 


84.41515 


-2.30742 


14.14 ± 





03 


13.93 ± 





03 


13.87 


± 





06 


13.74 


± 





08 




ill 


154 


84.43077 


-2.44282 


14.48 ± 





03 


14.45 ± 





04 


14.35 


± 





08 


14.15 


± 





09 




III 


174 


84.44174 


-2.45832 


13.71 ± 





03 


13.63 ± 





03 


13.61 


± 





05 


13.53 


± 





08 




ill 


188 


84.45503 


-2.36186 


14.14 ± 





03 


14.17 ± 





04 


13.94 


± 





06 


14.06 


± 





13 




111 


207 


84.46257 


-2.32070 


14.13 ± 





03 


14.18 ± 





03 


14.14 


± 





07 


14.20 


± 





11 




III 


221 


84.46809 


-2.44730 


14.24 ± 





03 


14.22 ± 





03 


14.29 


± 





07 


14.10 


± 





12 




III 


236 


84.47451 


-2.48824 


14.48 ± 





03 


14.46 ± 





04 


14.07 


± 





07 


14.13 


± 





13 




111 



Note. — Table [2] is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown 
here for guidance regarding its form and content. 

Column 9: III=non-excess, II=thick disk, I=class I candidate, EV= evolved disk, TD=transition disk candidate, 
DD=Debris disk candidate 

Classified as III in [03] 

2 Classified as EV in 3331 

3 Classified as II in SQ[2] 

4 Classified as I in P~2l 

5 Possible giant in the J-H versus H-K color-color diagram 

"Radial velocity consistent for stars in the older population llJeffries et al.ll2006h 
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Table 3. Median SEDs and Quartiles of optically thick disk stars 



Wavelength 


Median 


Upper 


Lower 


stars 


//m 


log , 

5 1.24F; 






used 


0.44 


-1.506 


-1.254 


-1.654 


22 


0.55 


-1.002 


-0.842 


-1.274 


51 


0.64 


-0.738 


-0.588 


-0.928 


60 


0.79 


-0.308 


-0.220 


-0.376 


52 


1.24 


0.0 


0.0 


0.0 


64 


1.66 


-0.037 


-0.005 


-0.065 


64 


2.16 


-0.209 


-0.125 


-0.245 


64 


3.6 


-0.563 


-0.400 


-0.645 


64 


4.5 


-0.733 


-0.580 


-0.836 


64 


5.8 


-0.910 


-0.748 


-0.978 


64 


8.0 


-0.962 


-0.804 


-1.076 


64 


24 


-1.238 


-1.029 


-1.430 


58 



Table 4. Median SEDs and Quartiles of non-excess stars 



Wavelength 


Median 


Upper 


Lower 


stars 


fim 


log 






used 


0.44 


-1.360 


-1.206 


-1.550 


60 


0.55 


-0.863 


-0.610 


-1.166 


116 


0.64 


-0.597 


-0.400 


-0.870 


128 


0.79 


-0.204 


-0.144 


-0.324 


97 


1.235 











137 


1.662 


-0.057 


-0.037 


-0.073 


137 


2.159 


-0.257 


-0.241 


-0.273 


137 


3.6 


-0.748 


-0.733 


-0.777 


137 


4.5 


-1.038 


-1.007 


-1.083 


137 


5.8 


-1.306 


-1.278 


-1.348 


137 


8.0 


-1.693 


-1.660 


-1.723 


137 


24 


-3.259 


-3.198 


-3.431 


42 1 



Including stars in the mass range of IMTTS and HAeBe 



Table 5. Disk fractions in the a Orionis cluster 



Range 


Thick Disks 

% 


Thick+Evolved Disks 

% 


HAeBe 


3.7±3.7 


14.8±7.4 


IMTTS 


19.6±5.9 


26.8±6.9 


TTS 


31.1±3.8 


36.3±4.1 


BDwarfs 


29.7±9.0 


33.3±9.7 


Total 


26.6±2.8 


33.9±3.1 
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Fig. 1. — False-color image of the a Orionis cluster. It is a 3-color composite of IRAC images, 
3.6 (blue), 4.5 (green), and 8.0 (red) /im. The field is centering at the star o Ori. 



SMQmDQ.Qs 39nkD0.0s 3SmO0.0£ 
Right ascension 



Fig. 2. — MIPS 24 pLm image of the a Orionis cluster. The solid diagonal lines indicate 



the I RAC field. The lower-right panel shows a map of dust infrared emission (ISchlegel et al. 



1998 ) illustrating th e location of the cluster in the OBlb sub-association (ring-like structure; 
Briceno et al.l 120051 ); the isocontours indicate levels of galactic extinction (Av=0.2, 0.4, 0.6 



& 0.8 mag). The IRAC (small box) and MIPS (large box) fields are shown in the inset. The 
positions of the young stellar cluster (<1 Myr) NGC 2024 and the Horsehead nebula are 
also indicated in this figure 
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Fig. 3. — [5.8]- [8.0] versus [3.6] color-magnitude diagram. Objects with large FWHM in the 
radial profile measured in all IRAC bands are likely to be extended sources (crosses). The 
density of objects with photospheric colors ([5.8]-[8.0]~0) decreases drastically for sources 
with [3.6]>14.5, most objects below this limit (solid line) are extended sources. We selected 
objects above the sol id line, where the contamination from extragalactic sources is expected 
to be less than 50% JFazio et al.l 120041 ). 
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12 3 

[5,8]-[8.0] 

Fig. 4. — IRAC color-color diagrams illustrating the color selection used to reject PAH-rich 
galaxies. Symbols are the same as in Figured Objects located redward from the solid lines 
are likely to be galaxies or faint companions with photometry affected by the primary stars 
of the stellar system. The open square is the star HD294268 (S0411) 
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Fig. 5. — 2MASS J magnitude distribution: the open histogram represents all the 2MASS 
sources located in the IRAC field. Striped and solid histograms represent the 2MASS sources 
detected in IRAC and MIPS, respectively. Our completeness limit for sources with [3.6]<14.5 
(see Figure [3} is J~14. Most objects with J<11 have MIPS detections. 
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Fig. 6. — V versus V-J color-magnitude diagram, illustrating the selection of members of 
the a Orionis cluster. Open circles represent stars in Sample 2. Solid symbols represent 
previous known members: photometric members (filled triangles) , spectroscopic members 
(filled squares), and members selected by iKharchenko et al.l (120041 ) combining photometric 
and astrometric data (filled circles). Stars with low membership probability combining pho- 
tometric and astrometric dat a are labeled with the p lus symbol. Open s quares represent t he 
X-ray sources in the sample (iFranciosini et al.ll2006l ). The ZAMS fr o m ([Siess et al.l 120001 ) is 
represented as solid line and the 3 Myr isochrone from lBaraffe et al.l (119981 ) is represented as 
a dashed line. The arrow represents the reddening vector (Av=l). The membership limit 
is defined by the dotted line; stars located blueward from this limit are rejected as likely 
non-members of the o Orionis cluster [See electronic edition of the journal for a color version 
of this Figure] 
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Fig. 7. — J versus J-K color-magnitude diagram, illustrating the selection of members of the 
a Orionis cluster. Symbols are similar to Figure [61 Stars rejected as likely non-members in 
Figure [6] are represented with crosses. The arrow represents the reddening vector (Av=l). 
Contamination by non-members are expected for stars with J>13. We reject as non-members 
of the a Orionis cluster stars located blueward from the dotted line. Stars selected using 
only the criterion defined in this plot are cataloged as uncertain members [See electronic 
edition of the journal for a color version of this Figure] 
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Fig. 8. — IRAC color color diagrams, [5.8]-[8.0] versus [3.6]-[4.5] (left panels) and [5.8]-[8.0] 
versus [4.5]- [5.8] (right panels). Members are plotted in the upper panels and uncertain 
members are plotted in the lower panels. Da shed box represents the loci of CTTS with 
different accretion rates (jD'Alessio et al.ll2005bl ). Open circles, filled triangles, filled squares 
and filled diamonds represent IRAC class III stars, stars with weak IRAC excess, IRAC class 
II stars and IRAC class I candidate (classified in §4.21) . respectively. Objects surrounded by 
open squares are transition disk candidates (see Figure [9]) [See electronic edition of the 
journal for a color version of this Figure] 
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Fig. 9. — [J-H] versus K-[24] color- magnitude diagrams for members (upper panel) and 
uncertain members (lower panel) in the a Orionis cluster. Symbols are similar to Figure [SJ 
Dotted lines define the locus expected for non emission stars. The dashed line at K-[24]=3.3 
(the color of the debris disk star f3 Pic) indicates the limit between optically thick disk 
objects and evolved disk objects; transition disk candidates are stars classified as IRAC class 
III or evolved disk objects located redward from the dashed line (objects surrounded by open 
squares). In general, the IRAC class I candidates are locate redward from the optically thick 
disk objects. The arrow represents the reddening vector (Ay=2) [See electronic edition of 
the journal for a color version of this Figure] 
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Fig. 10. — IRAC SED slope versus [8.0] for members (filled symbols) and uncerta in members 
(open symbols) in the a Orionis cluster. Dotted lines represent the limits from lLada et al. 
(120061 ) for different objects: IRAC class III stars (circles), IRAC class II stars (squares) and 
evolved disk objects (triangles). The dashed line represents the a median for IRAC class 
III stars The solid lines represent 3 a values around this median; we adopt the solid upper 
line as the boundary between IRAC class III stars and evolved disk stars. The dotted line 
at a=0 indicates the boundary between IRAC class II stars and objects with flat or rising 
SED, class I candidates (diamonds). Class III stars and evolved disk objects surrounded by 
open squares are transition disk candidates, since they show strong 24 /im excess consistent 
with having an optically thick disk [See electronic edition of the journal for a color version 
of this Figure] 
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Fig. 11. — IRAC SED slope versus J magnitude for stars in the a Orionis cluster. Symbols 
are as in Figure [TQl Vertical d ashed lines ind i cate approximate mass ran ges defined using 3 
Myr evolutionary models from lBaraffe et al.l (119981 . for <1 M ) and from lSiess et al.l (120001 . 
for >1 M ), assuming a distance of 440 pc and a reddening of E(B-V)=0.05 [See electronic 
edition of the journal for a color version of this Figure] 
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Fig. 12. — Disk fraction versus J magnitude. Each point represents a different mass range 
(see Figure fTTT) : from left-to-right, HAeBe, IMTTS, TTS and Brown Dwarfs. Solid circles 
+ dotted line indicate thick disks fraction. Open circles + dashed line are the fraction of 
disk-bearing stars. The lowest disk fraction is observed in the HAeBe mass range, while the 
highest frequency of disks is observed in the TTS mass range. 
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Fig. 13. — Spectral energy distributions for selected stars, illustrating the diversity of disks 
found in the a Orionis Cluster. The first 5 rows show stars in the mass range of TTS. Dotted 
and dashed lines represent the median SED of optically thick disks and non-excess stars of 
the a Orionis cluster, respectively. Error bars denote the quartiles of the optically thick 
disk median SED. Dot-dashed line, shown in the first panel, represents the median SED of 
optically thick disks in Taurus. The last row shows stars in the mass range of HAeBe; these 
are compared with the pho tosphere of a star with a similar spectral type (long-dashed line; 
Kenyon fc Hartmannlll995l ). Each panel is labeled with the internal running identification 
number and the classification from §4.21 
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Fig. 14.— Fraction of stars with near-infrared disk emission as a function of the age 
of the stellar group. Open circles represent the disk frequency for stars in the TTS 
mass range, d erived using JHKL observations: N GC2024, Trapezium, NG C2264 and 
NGC2362 from lHaisch et all (|200ll ), Chamaleon I fromlGomez fc Kenvonl fepQlh . and OB la 



and OBlb sub-associations from Hernandez et al.l (120051 ). Solid symbols represent the 
disk frequency calculat ed for stars in the TTS mass ran ge using IRAC data: Ta urus 



jHartmann e t alJ l2005ch. NGC7129 jGutermuth et ali l2004h. IC 348 Lada et al. 



37 and NGC7 160 (ISicilia-Aguilar et al.ll2006l ). Upper Scorpius ( jCarpenter et al. 
Chameleontis ( Megeath et al. 2005al ) and a Orionis cluster (this work) 



2006), Tr 



20061 ). T] 
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Fig. 15. — IRAC SED slope distribution for young stellar groups with differing ages. The 



dotted line defines the limit between thick disks and evolved disks. Fractions of disk- 



stars a nd stellar ages a re been included in each panel: Taurus ([Hartmann et al 



jearm g 



2005d), 



IC348 JLada et dJhood ). a Orionis (this work) and Tr 37 Jsicilia-Aguilar et al.ll2006h . The 



median of a slope and its quartiles for disk-bearing stars in each stellar group are indicated 
with arrows and error bars, respectively. 
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Fig. 16. — Diagrams illustrating the space distribution of members, disk fractions and the 
density of the cluster. The lower-left panel shows the location of members in the a Orionis 
cluster. Symbols are similar to Figure El Dotted circles define radial distances (at 0.1, 0.2 
and 0.3 degree) from the center of the cluster and solid lines define the semi-circles East- 
West (vertical) and North-South (horizontal) used to calculate the fraction of disks in bins at 
several radial distance from the star a Ori. Upper left panel shows the fraction of thick disk 
stars (filled circles + dotted line) and disk-bearing stars (open circles + dashed line) in radial 
distance bins defined by the semicircles East- West. Similarly, lower-right panel shows the 
fraction of disks in radial distance bins defined by the semicircles North-South. Upper-right 
panel shows the density of members in radial distance bins defined by the semicircles East- 
West (filled squares+solid line) and by the semicircles North-South (open triangles+dotted 
line) [See electronic edition of the journal for a color version of this Figure] 
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Fig. 17. — Cumulative distribution of non-excess stars and optically thick disk stars as a 
function of X-ray luminosity. A Kolmogorov-Smirnov (KS) test shows a significance level of 
64% when comparing both populations. This suggests that the X-ray emission has a similar 
origin in non-excess stars (like WTTS) and stars with optically thick disks (like CTTS) 



